Effect of spin-orbit interaction on a magnetic impurity in the vicinity of a surface 
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We propose a new mechanism for surface-induced magnetic anisotropy to explain the thickness- 
dependence of the Kondo resistivity of thin films of dilute magnetic alloys. The surface anisotropy 
energy, generated by spin-orbit coupling on the magnetic impurity itself, is an oscillating function 
of the distance d from the surface and decays as 1/d 2 . Numerical estimates based on simple models 
suggest that this mechanism, unlike its alternatives, gives rise to an effect of the desired order of 
magnitude. 
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The original observation that the amplitude of the 
Kondo resistivity in thin films of dilute magnetic alloys 
depends on the thickness L of the film, llL|3] has at- 
tracted considerable attention 0, 0, IE So far 
the most promising explanation of this phenomenon was 
given by Ujsaghy et al. |E : They suggested that a 
magnetic impurity, such as Fe in Au, near the surface 
of the host metal is subject to a magnetic anisotropy 
that blocks the spin-dynamics responsible for the Kondo 
scattering within a given distance, L c , to the surface, 
if the anisotropy is larger than the Kondo temperature 
Tr- of the magnetic impurities. As a consequence, the 
Kondo resistivity Rk{L) of a thin film relative to that 
of a 'bulk-like' thick film, R}^ lk (L), can be estimated for 
L ^> L c as 
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where the factor 2 on the rhs accounts for the two sur- 
faces of the film. Fitting the experiments with this for- 
mula yields L c ~ 180 A in case of Au(Fe) thin films 0. 

Ujsaghy et al. 0, Q based their arguments on cal- 
culations for a 5/2 Kondo impurity embedded in a semi- 
infinite host with spin-orbit interaction on the host atoms 
only. They found a surface induced anisotropy energy, 

H anis = K(d)S 2 z , 

where S z is the z-component of the impurity spin oper- 
ator and the anisotropy constant K(d) is a function of 
the distance d of the impurity from the surface. Using 
lowest order perturbation theory in both the spin-orbit 
coupling constant A and the effects of the surface they 
derived an expression for K{d) which for large d decayed 
as 1/d. Moreover, their quantitative estimates of K(d) 
gave sufficiently large values to explain the experimental 
facts in their more detailed calculations ISj. 



Subsequently, Szunyogh and Gyorffy studied the 
problem using a material specific, parameter-free first- 
principles approach, namely, spin-polarized relativistic 
calculations based on the Local Density Approximation 
(LDA). In these calculations, the semi-infinite host was 
taken into account without significant approximations. 
Consequently, the spin-orbit coupling on the host atoms 
and electrons' scattering between the impurity and the 
host atoms were treated on equal footing and to all 
orders in the coupling strength. They found K(d) ks 
cos(2k' F d) / d 2 , where 2k' F is the length of a spanning vec- 
tor of the host's Fermi Surface. Moreover, the size of the 
anisotropy energy turned out to be too small by orders 
of magnitude. Clearly, the LDA results would eliminate 
the whole idea of explaining the size dependent Kondo 
effect by 'surface induced magnetic anisotropy', if they 
didn't suffer from the well-known weakness of the LDA in 
describing spin-fluctuations on the impurity. LDA calcu- 
lations also underestimate 'Hunds rule correlations' sys- 
tematically. 

In this Letter we present calculations which include dy- 
namical spin- fluctuations at the same level as Refs. [a,LD|j 
but go beyond their approach in that here we treat 
the spin-orbit coupling in the semi-infinite host non- 
perturbatively, and we also incorporate the effects of 
spin-orbit coupling on the impurity. We find that, while 
the host-induced anisotropy, as proposed in Ref. 0, is 
negligible, an improved treatment of correlations and the 
strong (typically ~ 1 eV) spin-orbit coupling on the mag- 
netic impurity can lead to a dramatic enhancement of 
the surface induced magnetic anisotropy for impurities 
with partially filled d-shell. As it turns out it is large 
enough to explain the experiments. Furthermore, this 
anisotropy has a simple physical origin: for partially filled 
d-shells each spin-state has also an orbital structure. In 
a given spin-state, electrons on the deep d-levels lower 
their energy by hybridizing with the conduction electrons 



through virtual fluctuations. However, in the vicinity of 
a surface Friedel oscillations appear and, therefore, the 
density of states available for the d-electrons to hybridize 
with depends on the orbital state, and, thus, on the spin 
of the impurity (see Fig.^). This mechanism gives rise to 
an anisotropy that already appears to first order of the 
exchange coupling J between the magnetic impurity and 
the conduction electrons, decays as 1/cP and is orders of 
magnitude larger than the similar anisotropy induced by 
spin-orbit coupling on the host sites [|J. 

We shall first analyze the simplest possible model that 
captures all important aspects of the problem, and start 
with a magnetic impurity in a d 1 configuration such as 
a V A+ or Ti 3+ ion embedded into a simple cubic (sc) 
lattice. In this case, by Hund's third rule, strong local 
spin-orbit coupling will lead to a J = 3/2 multiplet that 
is separated from the J = 5/2 state typically by an en- 
ergy of the order of ~ 1 eV. The advantage of this model 
is that under a cubic crystal field the J = 3/2 multiplet 
remains degenerate (Tg double representation), and no 
anisotropy is generated when the surface is absent. For 
the sake of simplicity, we assume that the host atoms 
form a (001) surface of a simple cubic lattice. We also 
assume that the conduction band of the host is domi- 
nated by s-electrons, which we describe in terms of a 
single-band nearest-neighbor tight-binding model. 

The impurity's J = 3/2 multiplet can only hybridize 
with those linear combinations of s-states on the near- 
est neighbor atoms which transform according to the Tg 
representation. Within the full 12-dimensional subspace, 
spanned by orbitals on six neighbors and two spin-states 
per site, two such (four-dimensional) sets can be con- 
structed. One of these has p-type orbital structure and, 
therefore, hybridizes only weakly with the impurity's d- 
level. The d-type set, s m (m = — 3/2, . . . , 3/2), which 
hybridizes strongly with the J = 3/2 multiplet, splits 
into combinations with two different orbital characters 
as depicted in Fig. ^ the combinations s±i/2 couple to 
the impurity states J z — to — ±1/2, while s± 3 / 2 couple 

to J z — TTL — ±3/2. 

Assuming that the impurity-host interaction is mainly 
dominated by quantum fluctuations to the (non- 
degenerate) dP state, in lowest order of the hybridization, 
a Coqblin-Schrieffer transformation leads to the follow- 
ing effective exchange interaction, 01 
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where §to) stand for the four states of the Tg impu- 
rity multiplet, an d s m are creation and annihilation 
operators acting on the host states, respectively, and J 
denotes the effective strength of the coupling. We then 
employed Abrikosov's pseudo-fermion representation [ll| 
to calculate the splitting of the four states up to second 
order in J |6(. The first and second order contributions 
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FIG. 1: Sketch of the orbital structures of the d-type Tg 
combinations formed by s-orbitals at the nearest-neighbor 
sites around an impurity in a simple cubic host. The s±3/2 
combinations are composed from orbitals located at neighbors 
within the same plane, whereas the s±i/2 combinations have 
significant out-of-plane contributions. 



to the self energy at T — are given by 
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respectively. Here p m m'(e) denotes the local spectral 
function of the host computed in the absence of the ex- 
change interaction, i.e., J = 0, and Ep is the Fermi 
energy. [12J To compute g mm ' (e) , we made use of the so- 
called surface Green's function matching procedure [l3| 
that completely accounts for the semi-infinite geometry 
of the host. Note that Eqs. © and QJ, in addition to in- 
corporating quantum fluctuations of the spin at the same 
level as Ref. gg, also take into account the semi-infinite 
nature of the host non-perturbatively through the spec- 
tral functions (lif. 

In this simple case, tetragonal symmetry of the sc(001) 
surface implies that g mm i (e) , consequently, S mm / are di- 
agonal in to, to'. From time reversal symmetry it further 
follows that the states || ±|) and || ±|) remain de- 
generate. Thus the fourfold degeneracy of the J = 3/2 
multiplet is split by an effective anisotropy term, 



= KJi 



(5) 



with K = (S 3/2 - Ei/ 2 )/2. 

A clear understanding of the level splitting due to the 
vicinity of a surface obviously emerges from Fig. ^ the 
spectral density related to S3/2 orbitals extending in a 
single atomic plane differs from that corresponding to 
S1/2 orbitals that in fact take an average over three ad- 
jacent atomic planes. This is demonstrated in Fig. [2 
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FIG. 2: Spectral density functions corresponding to orbitals 
m = ±| (solid line) and m = ±^ (dashes), see also Fig. 
at the fifth layer below the (001) surface of a sc single-band 
host metal. The lattice parameter was chosen to 2.6 A and 
the hopping to 1.4 eV. 



where §3 (e) and g i (e) are plotted at the fifth atomic 
plane counted from the surface. In our numerical imple- 
mentation a layer spacing of 2.6 A (that corresponds to 
the atomic volume of fee Cu), and a hopping parameter 
V = —1.4 eV were chosen. Evidently, 03/2 (e) oscillates 
strongly within the energy range — 2\V\ < e < 2\V\ (the 
center of the valence band is set to zero), while £1/2 (e) 
behaves smoothly. 

Departing from the surface, the oscillations of g 3 / 2 (e) 
get more and more rapid while decreasing in magnitude. 
For large d, £1/2 (e) approaches £3/2 (e) and both tend 
to g bulk (e). Interestingly, at a given energy, they also 
display Friedel oscillations, 0] 



Q m (e;d) ~ A m cos(2k z (e)d + <j>)/d , 



(6) 



where k z [e) is the length of the extremal wave-vector 
parallel to the surface normal of the constant-energy sur- 
face in reciprocal space and 4> = or ir/2. Performing 
the energy integrations in Eqs. and Q, both self- 
energy contributions yield oscillations with a period of 
Tr/k z (EF) for large d and an amplitude ~ 1/d 2 , as also 
follows from an asymptotic analysis analogous to Ref. 9. 
The impurity's first and second order level splittings, 
E 1/2 , are plotted in Fig.^Jfor E F = V = - 1.4 eV. 
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Here we used J = 1 eV, a typical exchange coupling for 
Kondo impurities with a Kondo temperature of the order 
of a few Kelvins. As can also be obtained from analytic 
calculations, the period of the oscillations is 3 atomic 
layers (7.8 A) in this case. Remarkably, the second order 
self-energy diagram contributes about the same amount 
to the level splitting as the first order one. 

Note that the present theory also predicts that there 
will be impurities with nearly vanishing level splittings. 
However, for incommensurate Friedel oscillations the dis- 
tribution of A = |S 3 / 2 (d) — S 1 / 2 (c?)| within an interval 
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FIG. 3: First order (upper panel) and second order (lower 
panel) contributions to the level splitting of a d 1 impurity 
with strong on-site spin-orbit interaction as a function of the 
distance d from a sc(001) surface. See for numerical parame- 
ters in the text. For both cases the solid lines depict functions, 
Asm(2ivd/do)/d 2 (do = 7.8 A), fitted to the calculated values. 



[d — 8d, d + Sd] of a few atomic layers is peaked around the 
maximum values of the anisotropy within that interval, 
A max (d). Therefore, for large values of A max (c?) Tk 
(i.e., in the vicinity of the surface) only a very small frac- 
tion of the impurity spins will experience an anisotropy- 
induced splitting, A < Tk- As can be read off Fig. 
beyond d ~ 100 A the amplitude A max is still in the 
range of few tenths of a meV (1 K ~ 0.09 meV), and the 
typical level splitting is close to the values needed to sup- 
press the Kondo effect in thin films of alloys with Tk ~ 
0.1 - 1 K such as Au(Fe). 

As was stressed in the introduction, the above mecha- 
nism is a combined consequence of 'spin-orbit coupling' 
and 'Hunds rule' correlations at the magnetic impu- 
rity. Since the latter does not occur on the host sites 
we expect a much reduced contribution to the total 
anisotropy energy from these. In fact we would expect 
that, though the 'impurity spin' is treated classically, the 
host-induced anisotropy is described reasonably well by 
the first-principles LDA calculations, hence, the results 
for the anisotropy energy in Ref. |9j can be regarded as 
of the correct order of magnitude. This suggests that in 
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FIG. 4: LDA magnetic anisotropy energies for an Fe impurity 
in Au host as a function of it's layer position measured from 
the (001) surface. Squares refer to the experimental lattice 
constant (7.68 A), while triangles to the case when this lattice 
constant has been uniformly expanded by 50 %. 



models, such as that of Ujsaghy et al., [E 0> 13 where 
the spin-orbit interaction is restricted to the host atoms, 
the 'surface induced magnetic anisotropy' should be neg- 
ligibly small. From this point of view the fact that they 
find a level splitting sufficiently large to explain the ex- 
perimental data ought to be taken probably as a result 
of their largely analytic approximations. In fact, we also 
studied a somewhat modified version of our model where 
only host-induced spin- anisotropy appears and found, as 
in the first first principles calculations 0, that the ob- 
tained anisotropy was orders of magnitude smaller than 
the one discussed above. 

To lend further credence to the new mechanism of sur- 
face induced anisotropy we have attempted to simulate 
the effects of Hunds rule correlations on the first princi- 
ples calculations in Ref. 3. In calculations based on model 
Hamiltonians one can manipulate the parameters of the 
theory, like the hopping integrals or the impurity-electron 
interactions parameter, </, to desirable ends. Evidently, 
there is no such freedom in first-principles calculations, 
however, approximate. Nevertheless, in order to inves- 
tigate correlations between different effects one can con- 
strain the outcome of such calculations in other ways. 
As a quite crude device, we increased the lattice spacing 
of the Au host to mimic the band narrowing effects of 
correlations not included in LDA. In Fig. 0] we compare 
the surface induced anisotropy energies of an Fe impurity 
as calculated for the experimental lattice spacing with 
those for a lattice constant artificially enhanced by 50 
%. Note that because of the different lattice constants 
this comparison is made in terms of the impurity's posi- 
tion measured in the layer index rather than in physical 
distances from the surface. As expected from the consid- 
erations above, increasing the lattice spacing resulted in 



a dramatic, order of magnitude increase of the predicted 
anisotropy. Reassuringly, this result is quite robust since 
it occurs in spite of the fact that, beside the band nar- 
rowing required to increase the effect of the on-site spin- 
orbit coupling, such an increase of the distance between 
the atoms should also reduce the coupling between the 
impurity and the conduction electrons (J) and, hence, 
the anisotropy energy. Clearly, the observed overall en- 
hancement supports our contention that the origin of the 
large anisotropy energy and long critical length, L C1 is, 
indeed, the electron-electron correlations which are cor- 
rectly captured by the basic ansatz of our d 1 model but 
are neglected in the first principles LDA calculations. 
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